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"^ . Abstract. We construct higher order rogue wave solutions for the Gerdjikov-Ivanov equation 

On \ expHcitly in term of determinant expression. Dynamics of both sohton and non-sohton solutions 

is discussed. A family of solutions with distinct structures are presented, which are new to the 

Gerdjikov-Ivanov equation. 
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1. Introduction 

^ ■ Rogue wave is one type of natural disasters first found in deep ocean wave. Later, people 

QQ . observed similar phenomena in other physical breaches, such as optical physics, plasmas, cap- 
lO ! illary waves and so on [l-5j. There is a consensus that these rogue waves are the result of 
modulation instability waves. Meanwhile, breather solution usually comes from the instability 
^ ■ of small amplitude perturbations that may grow in size to disastrous proportions. Therefore, in 
^^r^ . mathematical understanding, rogue wave can be treated as a limit case of Ma soliton when the 
space period tends to infinity , or of the Akhmediev breather as the time period approaches to 
infinity [6j. Due to both theoretical frame and reality application, it is imperative to do further 
k3i ■ study about rogue waves, especially, higher order rogue waves of different models. 
^ . The Peregrine soliton which is located in time-space plane is one of the formal ways to explain 

the surprise phenomenon mathematically [7]. The solution appears from a non-zero constant 
and disappears to the constant background as time approaching infinity, but it develops a 
localized hump with peak amplitude three times of average waves in the intermediate times. 
The object of higher order rogue wave is now under intense discussion by different approaches 
such as iterated Darboux transformation method [8l[9], the algebra-geometric means [TO], and 
generalized Darboux transformation [llj. Recently, the expression of the first order rogue wave 
solution and the figure of the second order rouge wave for the Gerdjikov-Ivanov (GI) equation 
was provided by the third and fifth authors of present paper [12]. However, the higher order 
rogue wave solutions for this equation have not been studied. The main aim of this paper is to 
discuss higher order rogue wave solutions and describe their different structures. 

The nonlinear Schrodinger equation is one of the most important equations in physics, which 
can be derived from Ablowitz-Kaup-Newell-Segur system p!5l[T6 j. Considering the higher order 
nonlinear effects, the derivative nonlinear schrodignger equation with a polynomial spectral 
problem of arbitrary order [17j is regard as a model in a wide variety of fields such as weakly 
nonlinear dispersive water waves [TSj , nonlinear optics fibers [T9H2T] , quantum field theory [^22] , 
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plasmas [23] • The DNLS equations have three generic deformations, the DNLSI equation [24j : 

i<?t - Qxx + i{q^q*)x = 0, (1) 

the DNLSII equation [25j: 

kt + qxx + m*Qx = 0, (2) 

and the DNLSIII equation or the GI equation [26j: 

1 



iqt + qxx - ^q^ql + 2^^1*^ 



0. 



(3) 



In many circumstances, Darboux transformation has been proved to be one powerful methods 
to obtain sohton solutions [T3^jT!4^, breather solutions and rational solutions. The Darboux 
transformation and its determinant expression for the GI equations have been given in p!2l[27] . 
But there is not a straight extension construct higher order rogue wave solutions at the same 
eigenvalues. In this paper, we take a limit technique with respect to degenerate eigenvalues 
and Taylor expansion in Darboux transformation [lll[28[[29]. Based on this explicit method, 
we can further discuss the structure of solutions, both soliton and rogue wave. 

This paper is organized as following: In section 2, we review the general n-fold Darboux 
transformation for the GI equation. In section 3, explicit solutions are constructed, such as 
soltion, breather, position solutions and higher order rogue wave with two parameters Di 
and D2. By choosing different values of Di and i^2, we show four basic models, fundamental 
pattern, triangular structure, modified triangular structure and ring structure, and display their 
dynamical evolutions respectively in section 4. The conclusions and discussions are contained 
in the final section. 



2. Darboux transformation for the Gerdjikov-Ivanov equation 

In this section, we start with Lax pair of ([3]) to construct Darboux transformation. Consid- 
ering the spectral problem 

d^ij = {JX^ + QiX + Qo)iJ = U^j, 

^ dt^ = (2 JA^ + VsX^ + V2X^ + ^lA + Vo)^ = V^. 

where 



(4) 



xP = 



V 



4>{x,t, A) 
(plx,t,X) J ' 



J = 



0' i ' ^^ - ( r • ' ^° - 



-iqr * 

' 1. 



V 







-tqr 



J 



Vs = 2QuV2 = Jqr,V,= 



iq^ 
-ir^ 



Vo = 



-{rq,: - qr^) + -iq^r^ 







V 







-^(rqx-qrx) - -iq^r^ 



here A G C, ^/^ is the eigenfunction of (jlj) corresponding to the eigenvalue A. By the condition 

Ut-V, + [U, y] = , we get 

1 



igt + qxx + iq^r^ + -gV^ = 0, 



1 



(5) 



in - r^ 



ir Qx q r = 



This system admits the reduction r = —q*, and ([5]) becomes just one equation which is the GI 
equation ([3]). 
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From gauge transformation, we can construct new solution from initial data, i.e, if there exist 
some non-singular T, such that 



yW = {Tt + TV)T-\ 



(6) 



where t/'^l and y'^1 have the same form as U and V with q and r replaced by certain g'^1 and 
r'^l. Therefore, it is crucial to find an algebraic formula for T instead of the ([6]). 



2.1 N-fold Darboux transformation for the GI system 

The n — fold Darboux transformation for q^^\ r'^' of the GI system has been given in in [T2]. 
Since we need this result for our paper, we cite the main theorem as follows. 



Theorem 1. Let ^^ = 

problem , then (q^^\r^' 



Here, (1) for n = 2k, 



1,2, 



, n) be distinct solutions related to Xi of the spectral 



given by the following formulae are new solutions of the GI system. 



q + 2i— — , r^ ^ = r — 2i— — . 



^ 



11 



A29^2 



0n K^n 



Q 



12 



n 



21 



Q 



22 



(2) forn = 2k + 1, 



n^ = 



Q 



12 



Q 



22 



n-2^ 



Ai (fi Xi 

^2 9^2 ^2 < 



-A?(/)i 
-A^(/)2 



n-2. 



ArVn Ar'0n -a;:^^ 
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.. Ar 
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ArVi 
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A202 ■ 


■ Ar 
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Ar 


_9 

V2 


ArV2 


^n 


An0n • 


'• Ar 
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'Vn 


ArVn 


^1 


Al01 . 


.. Ar 


-'01 
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-Vi 
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A202 • 
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-V2 


-A^02 


^n 


An0n • 


• • Ar 


-'0n 
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-A::0n 
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(8) 



(9) 



n 



12 



Q 



21 



Q 



22 



(fi Ai0i 

^2 A2(/>2 

01 Ai^^i 

02 A2(/:^2 

01 Ai(/Pi 

02 A2(/:'2 

0n A^99^ 



\ 01 Ai 9^1 

\n—3i \n—2,^ 

X2 02 A2 (^2 



¥^2 



Ar^ 



Vn— 3 , 



\n-2.^ \n-l.^ 

^n ^n K ^' 



\n— 3,^ \n- 

Ai (^l \i 



\n— 3,^ 



ArV: 



-A^<^2 



ArVn Ar^0n -a;:^. 



Ar'^1 Ar'0i ArVi 

ArV2 Ar'02 ArV2 



ArVn Ar20^ Arv^ 



To solve the GI equation, we need to preserve the reduction condition r = — g*, i.e, under 
n steps of Darboux transformation, r'^l = — g'^l . Therefore, we can choose eigenfunctions 

0fc 



^k 



as foUows [I2J: 



(1) \k = -XI, and ^pk = 

(2) X2k = -As/e-i, and ^2k 



^k ) \^l 

^2k\ ^ f -^2k-l 
^2k J V ^2fc-l 

By this choice, we can construct sohton, breather, position solutions of the GI equation. But 
to get higher order rogue wave, we need to modify above Darboux transformation. At that 
time, the eigenfunctions and eigenvalues are no longer arbitrary, and only condition (2) works. 
We will discuss them in detail in the next section. 



3. Solutions of the Gerdjikov-Ivanov equation 
3.1 Solution with vanishing boundary condition 

We start with the soliton solutions by ([7j). Let g = 0, solutions of the spectral problem (|3j) 
with eigenvalues Xk are solved as 

^fc = ( M , 0fc = exp{-i{Xlx + 2Xtt)), ipk = exp(i(A2x + 2Xlt)). (10) 

Example: 

• n = 1. A^ = -Ai,Ai = i/3i and t/ji ^ ' 



1 



then 



czW=2Aexp-2'/'i'(--+2/3i2t) 



This is a plane wave with constant amplitude. 
n = 2. 
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1 K/, _ f -^ 



(1) Let A2 = AT, Ai = «! + i/9i and ipi = { ] fp2 = { T*^ ) 1 then 

\ "Pi J V v>i / 

gPl = -8 ""^^^^ (11) 

with 

Fi = - 2 i (tti^x + 2 aiH - 12 tti^/^i^t - ^i^x + 2 /5/t) , 

F2 =4tti/3i (4^1^ - 4t/3i2 + x) . 

It is a Hne sohton, and its trajectory is 

X = -At {ai - /5i) («i + /5i) , 

on the (x, t) plane. Let cti — )■ in (ITTl) . we obtain a rational travelling solution 



,[2] = zL ^^^ 

1- 41/31^0; + 16 i/3i^t' 



g'^'=4 /^V^,\^^.;,, (13) 



with an arbitrary real constant /3i. Its trajectory is defined explicitly by 

on the (x, t) plane. The above two solutions is plotted in Fig. [l](a) and Fig. [I](b). 
(2) Let Ai = — A^ = i/3i, A2 = — A2 = 1/^2 , then ipi = i ^l ) , ^^2 = ( ^l ) . we obtain 
a soliton solution 

^[2] = _2 ' [ P2 ^Pl ) ,^^. 

g-i(/5i-/52)(/3i+/32)(2t/3i2-x+2/32't)n _Q^i{(3i-(32){(3i^(32){2t(3i^-x^2(32h) ^ ^ 

Its trajectory is defined explicitly as 

in the (x, t) plane. It is plotted in Fig. [11(c). 

• n = 3. Let A^ = —A*, i = 1, 2, 3, we get a one-soliton solution. While a Ai = — A^, 
A3 = A2, we obtain a quasi-periodic solution. Since the analytic expressions are clumsy 
and tedious, we plot a particular quasi-periodic solution in Fig. 2(a). 

• n = 4. Under the different reduction conditions, we obtain three kinds of solutions. One 
of them is a soliton solution under the quasi-periodic background, the others are second 
soliton solutions. A soliton solution under the quasi-periodic background is plotted in 
Fig. 2(b). 

Remark. If the eigenvalues share the same value, T becomes degenerate. Next, we will 
consider the degenerate cases by limit technique. 

Given ip = '0(A) , define ip[i^j^ k] by 

X^ij = il,[i,j,0]+ij[i,j,l]e + ij[i,j,2]e' + --- + i,[i,j,k]e'' + --- , (15) 

with 

^K J, k] = ;^^((A^ + eYiPiX, + e)). 
In other words, '0[i, j, k] is the coefficient of e^ if we expand X^ip at A = A^ by Taylor expansion. 



Now, we consider the degenerate case of n = 4. According to the formula, 

2i5n 



7[41 = 



5 



12 



(16) 



with 



Su 



(f)i Xixpi Ai^(^i -A/(^i 

^2 X2tp2 A2^^2 -A2''^2 

4>3 hi'3 h'^4>3 -h'^4>3 

(j)4 X41P4 A4 (j)4 — A4 (j)4 



(j)l Ai-^i AiVl Ai^t/^i 

02 A2'?/'2 A2^02 A2^'02 
-"12 = 2 3 

(^^is Ast/^s A3 (f)3 A3 -i/;3 

<^4 A4'04 A4 (^4 A4 'i/^4 

Notice that these eigenfunction vectors t/^j(« = 1,3) have the same formula except only 
different eigenvalue Aj. Thus, if we set Ai = A3, the determinants are linearly dependent. Then, 
we will obtain a trivial solution q^*^ = 0. To get non-trivial solution, we need to adopt the 
Taylor expansion technique. 

The main process: i) For the first(second) row, we can substitute the eigenvalue Ai(A2 = XI) 
and eigenfunction ipi{ilJ2) directly, ii) A3 — )• Ai(A4 — > A2), we expand the elements of the 
third(forth) row at A3 = Ai + e(A4 = A2 + e) by Taylor expansion at first, then subtract the 
first(second) row from the third(forth) row. iii) Taking e ^ 0, the terms with higher order of 
e will vanish. and the new expression of gW is obtained as the same form as eq. (I16|l . but with 
the values for ^n and 5i2 given by 



5n = 



0[1,O,O] (^[1,1,0] 0[1,2,O] 

(/)[2,0,0] ¥^[2,1,0] 0[2,2,O] 

(/)[!, 0,1] vp[l,l,l] 0[1,2,1] 

(/>[2,0,1] vp[2,l,l] 0[2,2,1] 



-0[1,4,O] 
-0[2,4,O] 
-0[1,4,1] 



,4,1] 



^2 



0,0] ¥^[1,1,0] 

0[2,O,O] ¥^[2,1,0] 

0[1,O,1] vp[l,l,l] 

(/>[2,0,1] ip[2,l,l] 



2,0] ¥^[1,3,0] 

2,0] ¥^[2,3,0] 

2,1] (/p[l,3,l] 

2,1] (/p[2,3,l] 



(17) 



Substituting Xi = ai + i/3i and eigenfunctions of eq. (ITOl) into the above formula, we obtain 
a positon solution 



[4] _ Gi + i(?2 
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(18) 



where 

G2 = - 16cos(F2)tt?(-32/3?tt?^sinh(Fi) + 16/3itt^tsinh(Fi) + 4/3?a;sinh(Fi) + 4tt^/3^a;sinh(Fi) 

- 48^ftsinh(Fi) - tti cosh(Fi))/5i, 
Gs =16 ai^/^i^a; + 64 cti'^t/Ji^ + 2 ai/^i^ sinh (2 F2) - 384 ai^/3i^t - 16 cti^/^i^x 

+ 2 ai^/?! sinh (2 F2) + 64 ai^l3iH, 
Ga = - 256 ai^x(3ih + 32 tti^a;^/?!^ + ai^ + 64 /^/o^^tti^ - 256 ai'^x^i^t + 256 tti Vi^^;^ 

+ 3072 (3ih'^ai^ + 2048 ^iH'^ai^ + 32 tt/a;^/?!^ + 512 ai^^t^(3i^ + 512 /5/°t2«i2 

+ 2048 aiH'^^i^ + /3/ + 256 (3i^xaiH + ai^ cosh (2 F2) - /5i^ cosh (2 F2) , 
Fi =4 «i/5i (a; + 4 toi^ - 4 t/^i^) , 
F2 = - 24 fti^t/^i^ + 4 /5/t + 4 aiH - 2 /5i2a; + 2 ai^o;. 

Its dynamical evolution is plotted in Fig. 3(a). 

If we let cti — > in above procedure, we will get the second order rational solution. With 
these parameters, the general solution can be given as following: 

^"' = n ^V exp(-2i/3i^ {-x + 21^^)), (19) 

with 

Gi = - ^Pi{-l2iPlx + IQ^i^lxH + AmQi^lH^ - QAi^lx^ - 48i/3i^t - 3072i/3f xt^ - 3 + 2?>{)U^Pl 

-768x/3i^t + 48/3iV), 
G2 =3 - 768x/3/t + 4096/3i^^tx^ - 24576 /3i^W + 65536 /3i^^xt^ + 4608t^/3i^ - 65536 /3i^^t^ 

-256/3iV + 96/3iV, 
G3 = - 12288 /3i^°xt2 + 16384 /3i^¥ - 256/3i^x^ - 48/3i^x + 3072/3i^x^t + 576 /3/t. 

The equations (IT81) and (IT91) represent the interaction of sohtons and rational sohtons, re- 
spectively. A simplify analysis shows that they possess phase shift when t -^ ±oc. Which is 
different from general 2-soliton solutions. They are shown in Fig. 3. 

This process works for general n. 



Theorem 2. For n = 2k, '0=1 1 ^5 eigenfunction vector of the GI system, the expression 

of n-positon solutions in terms of determinant is obtained: 

^["1 = 2i^, (20) 
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where 



ni 



^12 



0,0 


(^[1,1,0] •• 


¥p[l,n-3,0' 


<^[l,n-2,0' 


-(^[l,n,0' 


0,0 


(/9[2,1,0] ■• 


¥p[2,n-3,0' 


<^[2,n-2,0' 


-<^[2,n,0 


0,1 


(^[1,1,1] ■• 


¥?[l,n-3, 1 


<^[l,n-2, 1 


— (/)[l,n, 1 


0,1] 


(^[2,1,1] •• 


■ ¥p[2,n-3,l" 


<^[2,n-2,l 


-<i-[2,n, 1 



0[1,O,A;-1] (^[l,l,fc-l] 
0[2,O,A:-1] ^[2,1,A:-1] 



(^[1,1,0] 
^[2,1,0] 
^[1,1,1] 
^[2,1,1] 



0,0] 
0[2,O,O] 
0[1,O,1] 
0[2,O,1] 



0[1,O,A:-1] ^[1,1,A:-1] 
0[2,O,A:-1] ^[2,1,A:-1] 



ip[l,n-?,,k-l] 0[l,n-2,A;-l] -<^[1, n, fc - 1] 
(/7[2, n - 3, A: - 1] 0[2,n - 2, A; - 1] -(/)[2, n, A: - 1] 



ip[l.,n — 3,0] 
(/7[2,n-3,0] 
(/7[l,n-3, 1] 
</p[2,n-3,l] 



<^[l,n 
(/)[2,n 
^[l,n 
(t>%n 



2,0] 
2,0] 
2,1] 
2,1] 



(^[l,n-l,0] 
^[2,n-l,0] 
^[l,n-l,l] 
^[2,n-l,l] 



(/7[l,n-3,A:- 1] (/)[l,n - 2, A; - 1] ^[l,n - 1, A: - 1] 
(/7[2, n - 3, A: - 1] 0[2, n - 2, A; - 1] ^[2,n - 1, A: - 1] 



Proof. For the entries in the first column of Hn ( 

01 = '/'[I, 0,0], 

02 = <^[2,O,O], 

03 = (/)[l,O,O] + (/)[l,O,l]e, 

04 = <^[2,O,O] + 0[2,O,l]e, 



(j)n-i = <^[1, 0, 0] + <^[1, 0, l]e + 0[1, 0, 2]e2 + • ■ • + 0[1, 0, A; - l]e'-\ 
<^„ = <^[2, 0, 0] + <^[2, 0, l]e + 0[2, 0, 2]e2 + • ■ • + 0[2, 0, A; - l]e'=-\ 

Taking the similar procedure to the other entries in r2ii,l]i2. FinaUy, the gt"' can be obtained 
through simple calculation and certain limit. D 

Let «i ^ in above formula, we will get the n-rational solution. 

3.2 Solutions with non-vanishing boundary condition 

In this section, we consider solutions from non-trivial seed, and this will give higher order 

rogue wave. In general, we start with q = ce''"^"*^" '^'^ «-V)*l, a,c & C^. Then the corresponding 
eigenfunctions ipk associated with A^ 



Mh)\^(DizuiiXk)[h k] + D2W2iXk)[l, k] - D2Wi*(A^)[2, A:] - Dit^*(A^)[2, A:] 
ifikiXk)] Uit^i(Afc)[2, A:] + D2ZU2iXk)[2, k] + D2t^i*(A^)[l, A:] + Diiu^(A^)[l, A;] 



(21) 



where 



zui{Xk)[2,k] 



(ia-ic^+2iAfc^+2ci) 
2AfcC 



exp( — Ci(ta — X — Z,LAk 

exp(— Ci(ta — X 



r)j.\ 2^ I iiax-{a'^+c^a- 



+ 



-)t). 



\ 



2t\k^) 



\{ax — {o?-\-c^a—^)t) ^ 



I 



W2{\k)[^,k] 



exp(ci(ta — X — 2tXk^) + 



i(ax — (a'^-\-c'^a—^T-)t) ^ 



^ (iQ-ic^+2iAfc^-2ci) 
2AfcC 



exp(ci(ta — X — 2tXk^) 



i(ax — (a'^-\-c'^a—^)t) n 



"^'^^'^ - [ miXk)[2,k] ) ' "^'^^'^ - [ ^2iXk)[2,k] 



Cl 



V-c^ + 2c2 



a — a^ 



4aAi 



4A, 



(22) 



In the following, we set L^i = 1 and D2 = I- Under the circumstance, Xu et. al. [12j have 
obtained two kinds of breather solutions. One is time periodic breather, and the other is space 
periodic breather. Besides, the first order and the second order rogue wave solutions have also 
been obtained by a certain limit from breather solution. However, the method is diflBcult to 
calculate the higher order rogue waves. What we do here is to take the limit in the determinant 
expression of solution for the GI equation directly, which is the same method as we do for the 
NLS equation [29j. 



Proposition 1. Let n = 2k, iJj 



^ 



is eigenfunction vector, assuming 



Ai = -V^ 



2a 



-ic 



and 



A. 



1. 



- Vc2 -2a + -ic. 



the the formula of the n-th order rogue solution is obtained 



(23) 



where 



(^11 = 



,0,1] ^[1,1,1] 

0[2,O,1] ^[2,1,1] 

(/.[1, 0,2] (^[1,1,2] 

0[2,O,2] ^[2,1,2] 

(/)[l,0,k] ^[1,1, k] 

0[2,O,A;] ip[2,l,k] 



ip[l,n-3,l] 
^[2,n-3,l] 
(^[l,n-3,2] 



l,n-2,l] 
,n-2,l] 
,n-2,2] 



-(/)[!, n,l] 

-(/)[2,n,l] 

-^[l,n,2] 



^[2,n-3,2] (p[2,n-2,2] -(p[2,n,2] 



(p[l,n — 3,k] (f)[l,n — 2,k] —(f)[l,n,k] 
'^[2,n-3,k] (f)[2,n-2,k] -(/)[2,n,k] 
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n2 



,0,1] vp[l,l,l] 

^[2,0,1] ¥^[2,1,1] 

0[1,O,2] vp[l,l,2] 

(/.[2,0,2] vp[2,l,2] 

0[1,O,A:] ^[1,1, A:] 

0[2,O,A:] (^[2,1, A:] 



(/?[!, n- 3,1] (/)[l,n-2, 1] ip[l,n-l,l] 

¥?[2,n-3, 1] <^[2,n-2, 1] ¥7[2,n-l,l] 

(/7[l,n-3,2] (j)[l,n-2,2] ip[l,n-l,2] 

(/7[2,n-3,2] (/)[2,n-2,2] (/7[2,n-l,2] 



<^[l,n — 3, /c] 
i^[2,n — 3, k] 



l,n — 2,k] <^[l,n — l,k] 
,n — 2,k] <^[2,n — 1, k] 



Proof. For the entries in the first column of i^n®, 

0i = 0[l,O,l]e, 

02 = 0[2,O,l]e, 

03 = 0[l,O,l]e + <^[l,O,2]e2, 

04 = 0[2,O,l]e + ^[2,O,2]e2, 



^n-l 



1, 0, l]e + (/)[1, 0, 2]e^ + ■■■ + 0[1, 0, A:]e^ 
, 0, l]e + (/)[2, 0, 2]e2 + ■ • ■ + 0[2, 0, A:]e^ 



Taking the similar procedure to the other entries in r2ii,r2i2. Finally, the gt"' can be obtained 
through simple calculation and certain limit. 



D 



As application, we give some explicit solutions next. 

• first order rogue wave 
For n = 2, we get the first order rogue wave. 



,[2] = 9l 



exp ( - i (2 aa; - 2 ta^ - 2 tc^a + tc^) j 



(24) 



with 



G2 



8 c^aH^ + 12 a^cH^ + 8 c^aHx - 8 cHax - 2 c^ax^ 
2ic^x + 2 c^x^ - QicH + 2 cH^, 
8 c^aH^ + 12 a^cH^ + 8 c^aHx - 8 cHax - 2 c^ax^ 
2ic^x + 2 c^x^ + 2icH + 2 c^t^ 



12iac^t-6ac^t^-3 
AicHa - 6 ac^i^ + 1 



Use this method, we also get similar result of [12]. Moreover, our method is general 
to product higher order rogue wave. Assuming x ^ oc, t ^ oc, then g'^1 -^ c^. The 
maximum amplitude of g'^1 occurs at origin and is equal to 9c^. A first order rogue wave 
with particular parameters is shown in Fig. 4. 

• Higher order Rogue Wave 

Generally, the expression of rogue wave becomes more complicated with increasing n. 
For convenience, we use numerical simulations to discuss the higher order rogue wave. 
In the following, set a = and c = 1. 
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When n = 4, we can obtain the second order rogue wave according to the formula 

^"' = f ^-P (^i^) ' (25) 

with 
Gi =45 + 8 1^ - 24 it^x - 144 it^x^ - 48 ifx^ - 72 itx^ + 576 ifx + 288 itx^ 

- 144 t^x - 504 t^x^ - 144 tx^ - 24 ix^ - 528 it^ + 414 it + 90 ix - 72 it^ 

+ 48 ix^ + 24 t^x^ + 24 t^x^ + 504 tx - 198 x^ - 486 1^ - 60 x^ - 60 1^ + 8 x^ 
G2 =9 + 48 t^x - 216 t^x^ + 48 tx^ + 24 it^ - 24 ix^ + 24 itx^ + 48 it^x^ + 198 it 

- 24it^x + 180t^ + 24 tV + 24 tV + 8t^ + 8x^ - 12x^ + 666 1^ + 90x2 

- 72 tx + 336 it^ - 48 ix^ - 54 ix + 288 it^x - 48 it^x^ 

Besides, we can also obtain the fc-th {k = 3,4,5,6,7) order rogue wave solution of 
the GI equation. Since their analysis expressions are too cumbersome, we omit them. 
Their dynamical evolutions are displayed in Fig. [51 From the figures, we find these 
local peaks all own a high amplitude on their center, and there are many small peaks 
locating around the central peak. Compared with the higher order rogue wave of the 
NLS equation, the central high peak of higher order rogue wave of the GI equation 
looks like higher than the same order of the NLS equation. However, through detailed 
analysis, we find the amplitude of the k-th order rogue wave solution of the GI equation 
and the NLS equation are both 2(/c + l)c (c is the boundary condition of the seed 
solution) [29j. 

4. The dynamics of Rogue Wave 

In above section, we obtain the fundamental pattern of the higher order rogue wave solution 
for the GI equation with Di = 1 and D2 = I in ( l2T]) . Actually, Di and D2 are arbitrary 
constant (or go to constant). In this section, we set Di and D2 as following: 

Di =exp(-ici(5o + Sie + 826^ + ^3^^ + • • • + ^^-le^"^)), 

D2 =exp(ici(5o + Sie + 826^ + Ssc^ + • • • + Sk-ic^'^)). 

Here aSi, 5^2, . . . , Sk-i G C. Although the terms with nonzero orders of e vanish in the e ^ 
limit, analysis and numerics prove that their coefficients ^^^(z = 1, 2, . . . , fc — 1) have a crucial ef- 
fect on the dynamics of higher order rogue wave. In the following, our main task is to talk about 
how do these parameters control these different spatial-temporal structures at the same order k. 

4.1. Solutions with one parameter 

• Sq: Fundamental pattern 

When S'i = (i 7^ 0), we obtain a trivial translation. A special case for the first order 
rogue wave is shown in Fig. [61 It just change the location of the rogue wave. Actually 
we can shift the location of rogue wave to arbitrary position by hanging the value of 
Sq. The case for the NLS equation had been given in [30j. Moveover, the k order rogue 
waves have n{n + l)/2 waves. Starting from large negative t, n small peaks, then a 
row of n — 1 larger peaks etc., the central high amplitude wave appears. The process is 
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reversed in positive t. The dynamical evolution can be observed for the NLS equation 

distinctly in f3T] . 

Si: Two triangular structure 

(1) Triangular structure 
Let all each coefficient Si = except Si. The higher order rogue wave solution 
of the GI equation split into a triangular structure. This structure of k-th order 
rogue wave contains k{k + l)/2 first order fundamental patterns, which make up k 
successive arrays, and these arrays possess /c, /c — 1, • • • ,1 peaks respectively. These 
structures of the k-th [k = 2, 5, 7) order rogue wave are shown in Fig. [71 From 
the Fig. [Z](c), we observe seven arrays, and each of them has 7,6,5,4,3,2,1 lo- 
cal maxima respectally. And the orientation of these triangular structures in the 
(x,t)-plane remain the same. Our result is a general case of f32] . 

(2) Modified triangular structure 
Actually, the outer triangular is independent of the inner triangular in above tri- 
angular structure. For instance, when A; = 5, the outer 12 local maxima and the 
inner 3 fundamental patterns both make up a triangular, but these tow triangles 
are irrelevant. If we alter the appearance of ( I2T]) and mixing coefficients Di and 
D2 as following: 

Mh)\ ^ fD,zu,{X,)[l, k] + DiW2{Xk)[l, k] - D2wl{Xl)[2, k] - D2w;{Xl)[2, k]\ 
M>^k)J \DiWi{Xk)[2,k] + D^W2{Xk)[2,k] + D2Vul{Xl)[l,k] + D2w;{Xl)[l,k]J ' ^^'^ 

with 

Di =exp(-ic?(5o + Sie + ^2^^ + ^3^^ + • • • + ^^-le^"^)), 

D2 =exp(ic?(5o + Sie + S^e^ + S^e^ + " - + Sk-ic'-')). 

Remark: the expression of (l27l) is different with (|2TI) . 

The inner three peaks form a second fundamental pattern inversely. It is a new 
structure which has never been obtained for the NLS equation. Its dynamical evo- 
lution is shown in Fig. [51 For k-th order rogue wave, we conjecture that there are 
3k — 3 first order rogue wave locating the outer triangular shell and an {k — 3)-th 
order rogue wave locating in the center. 

Sk-i'. Ring structure 

If Si = {i y^ k — 1) ^ we obtain a ring structure. A high maximum peak is surrounded 
with some local maxima. By simple analysis, we find the peaks locating on the outer 
shell are all first order rogue waves, and the number of them is 2/c — 1. Besides, the inner 
peak is a higher order rogue wave (except fc = 3), whose order is fc — 2 (when fc = 3, 
the inner peak is a first order rogue wave). For example, when fc = 7, there are 13 first 
order rogue wave locating on the outer shell, and a fifth order rogue waves locates in 
the center. Some of these structures are shown in Fig. [91 

(1) Sk-i^ Si: Ring-triangle 

If Si is also non-zero, the inner higher order rogue wave will become a triangular 
structure. For instance, when fc = 5, the inner second order rogue wave is split 
into a triangle with 6^4 = 5 x 10^, Si = 500. Besides, when fc = 6, 7, the sim- 
ilar structure is obtained. Therefore, we can conclude reasonably that the inner 
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higher order rogue wave is able to spht into triangular structure when Sk-i is big 
enough and Si ^ 0. We display some of these special models in Fig. [101 Notably, 
the orientation of the triangle in 7-th ring-triangle model is different with the others. 

(2) 5fc_i, 5fc_3: Multi-ring 

Being similar to the above case, the inner higher order rogue wave can be split 
into a ring structure too. Keeping S}^-\ ^ and setting Sks 7^ (fc > 3), we 
will obtain a ring structure. In this case, if the central peak is still a higher order 
rogue wave, we can continue splitting it into first order fundamental models. For 
example, when k = 6 with S^ = 8 x 10^^ Ss = 80000, we obtain a multi-ring 
structure with a seconde order rogue wave locating in the center. Its dynamical 
evolution is shown in Fig. [TTl Under this circumstance, assuming Si ^ 0, then 
the second order rogue wave is split into a triangular structure. Its evolution is 
displayed in Fig. [121 Naturally, these similar structures are also able to found for 
k = 7. Since the third order rogue wave possesses two structures i.e. triangular 
structure and ring structure, the inner structure of the seventh can be both triangle 
structure and ring model. Their evolution is shown in Fig. [131 

In general, there are k — 1 free parameters for k-th order rogue wave solution, which possesses 
2^ combinations. But there are only 4 basic models, and by choosing proper parameters, we 
also get dynamical evolutions of the 4 basic patterns. Therefore, we can get a hierarchy of 
higher order rogue wave of the GI equation with combination structure of above four basic 
models. 

5. Conclusion 

In this paper, we modified the generalized Darboux transformation to get explicit solutions 
for the GI equation. In the case of vanishing boundary condition, we obtain soliton, rational 
traveling soliton, breather type soliton, and soliton colliding with breather type soliton. The last 
two kinds of solutions are new for the GI equation. Moreover, we get the formula by applying 
the Taylor expansion and limit technique, when eigenvalues share the same value. Under 
the condition of non- vanishing boundary condition, we give the formula of the A^-rogue wave 
solution for the GI equation in proposition 1. As applications, we give the expressions of higher 
order rogue wave solutions and discuss their structures. Further more, we generate solutions 
with different structures by adjusting the free parameters Di and D2- In summary, there are 
four basic patterns for higher order rogue waves of the GI equation: fundamental pattern, 
triangular structure, modified structure and ring structure. By choosing proper parameter, we 
can get solutions with combination structure of basic ones. 

Our results clearly show the connection between shift parameters and structures of rogue wave 
due to the explicit formula. These types of solutions are new to the GI equation. Moreover, 
some basic models such as the fundamental pattern, ring structure and triangular structure 
have been found in higher rogue waves for other equations such as the NLS equation [9j, but 
the modified triangular structure is new to the GI equation. Meanwhile, we get a family of 
solutions with combination structures. All of these will help us in founding univeral properties 
for rogue waves and better understanding of "waves appear from nowhere and disappear without 
a trice" [6j . 
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(a) 



(b) 



(c) 



Figure 1. (Color online) The dynamics of |g[^l| . (a) Eq. ( ITTI) soliton solution 
with ai = 0.2, /3i = 0.4. (b) Eq. (IT31) rational travelling solution with /3i = 0.4. 
(c) Eq. dUD soliton solution with /3i = 0.45, ^2 = 0.55. 





(a) 



(b) 



Figure 2. (Color online) (a) The dynamics of soliton solution quasi-periodic 
solution on the (x, t) plane, (b) soliton solution under the quasi-periodic back- 
ground on the (x, t) plane. 







(a) 




(b) 



Figure 3. (Color online) The dynamics of 2-posoliton solution and 2-rational 
solution on the (x, t) plane, (a) Eq. ([IHD with ai = 0.5, /3i = 0.5. (b) Eq. ([I8j) 
with ai -^ 0. 
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Figure 4. (Color online) The dynamics of the first order rogue wave solution 
|g[^l I Eq. ( !24l) on the (x, t) plane with a = and c = 1. The maximum amplitude 
of Ig'^l I occurs at t = and x = and are equal to 9. 




c^MS 



(a) 








(c) 



Figure 5. (Color online) The dynamics of higher order rogue wave, (a) The 
second order rogue wave, (b) The fifth order rogue wave, (f) The seventh order 
rogue wave. 





(a) 



(b) 



Figure 6. (Color online) The fundamental pattern of first order rogue wave. 
(a) The first order rogue wave with Sq = 5^ the maximum amplitude occurs at 
i = — 5 and x = 0. (b) The first order rogue wave with Sq = —5^ the maximum 
amplitude occurs at t = 5 and x = 0. 
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(a) 



(b) 



(c) 



Figure 7. (Color online) The triangle structure of higher order rogue wave, (a) 
The second order rogue wave with Si = 1000. (b) The fifth order rogue wave 
with Si = 800. (c) The seventh order rogue wave with Si = 250. 





(a) 



(b) 



Figure 8. (Color online) The modified triangular structure with a higher order 
rogue wave locating in the center, (a) An overall profile with Si = 1000. (b)A 
local central profile of the right panel. 






(a) 



(b) 



(c) 



Figure 9. (Color online) The ring structure of higher order rogue wave, (a) The 
third order rogue wave with 5^2 = 50000. (b) The fifth order rogue wave with 
5^4 = 5000000. (c) The seventh order rogue wave with Sq = 5 x 10^. 
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(a) 



(b) 
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(c) 



(d) 



Figure 10. (Color online) The ring-triangle structure of higher order rogue wave. 
(a) The forth order rogue wave with 5^3 = 5 x 10^, Si = 500. (b) The fifth order 
rogue wave with 5^4 = 5 x 10^, Si = 500. (c)The sixth order rogue wave with 
Sq = 5 X 10-^^, Si = 1000. (d) The seventh order rogue wave with 6^7 = 1 x 10-^-^, 
Si = 150. 




Figure 11. (Color online)The multi-ring model of the sixth order rogue wave 
with S's = 8 X 10^, 5^3 = 80000, the inner peak is still a higher order wave 
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(a) The entire structure 



(b) The inner structure 



Figure 12. (Color online) The multi-ring model of the sixth order rogue wave, 
the inner higher rogue wave is split into a triangle. 





(a) The entire structure 



(b) The inner structure 





(c) The entire structure 



(d) The inner structure 



Figure 13. (Color online) The two kinds of multi-ring structure of 7-order rogue 
wave solution with three parameters, (a) The 7-order rogue wave with Sq = 
5 X W^\ 5^4 = 1 X 10^ Si = 100. (b)The 7-th order rogue wave with ^e = 5 x W^\ 
S^ = lx 10^ ^2 = 1000. 
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